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We report experimental results on the dynamics of heavy particles of the size of the 
Kolmogorov-scale in a fully developed turbulent flow. The mixed Eulerian structure func- 
tion of two-particle velocity and acceleration difference vectors (<5 r v • <5 r a p ) was observed 
to increase significantly with particle inertia for identical flow conditions. We show that 
this increase is related to a preferential alignment between these dynamical quantities. 
With increasing particle density the probability for those two vectors to be collinear 
was observed to grow. We show that these results are consistent with the preferential 
sampling of strain-dominated regions by inertial particles. 



1. Introduction 

Turbulence occurs whenever fluid viscous forces are small compared to the dominant 
driving forces of the flow. In practice this includes most macroscopic natural and tech- 
nological flows. Very frequently, these turbulent flows are loaded with small, passive 
particles with non-negligible inertia, i.e. particles do not exactly follow fluid motion. 
Familiar examples of such flows include clouds in the atmosphere, sand storms in arid 
regions, flow driven transports of grains in pharmaceutical and food industries, fluidized 
beds in chemical engineering, and fuel sprays in combustion engines and turbines. The 
understanding of particle-turbulence interaction is therefore of central importance to 
our ability to make advances in key economical and societal issues like energy genera- 
tion, climate change, and pollution control. Comparing to the fluid turbulence problem 
itself, which is well described by the non-linear Navier-Stokes equations, the particle- 
turbulence interaction problem is on much weaker foundations. In principle, the particle 
dynamics can be obtained by solving particle equations coupled with the Navier-Stokes 
equations for the fluid with no-slip boundary conditions on particle surfaces. However, 



this approach is a formidable task for theoretical and numerical investigations (Happel 



& Brenner 1965 Homann & Bee 20101. Most theoretical analyses and numerical sim- 
ulations therefore must rely on simplified equations, in which the particle-turbulence 
interactions are modeled. An interesting result obtained from these simplified equations 
is that in a turbulent flow particles with small to moderate inertia and sizes smaller than 
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PpIps 


dp (fim) 


dp/v 


St 


Polystyrene 


1.06 


74 ± 10 


0.98 


0.08 


Glass 


4 


75 ±8 


1.04 


0.24 


Steel 


7.8 


75 ± 15 


1.04 


0.44 



Table 1 . Characteristics of the different particles used in the experiments. The Stokes number is 
defined as St = dp/12/3n 2 , where d p is the particle diameter, /3 = 3p//(2p p + pf) is the modified 
density ratio, that takes into account the added mass effect, and n is the Kolmogorov length 
scale of the turbulence. 



the Kolomogrov length scale tend to move to regions with high strain and low vortic- 
ity ( Maxey fe Riley|1983 ). This so-called "preferential concentration" of inertial particles 
is under intensive investigation (Falkovich fe Pumir||2004 Zaichik & Alipchenkov 2009 



Chun et «Z.|[2005| |Collins fc Keswani||2004| |Sundaram fc Collins||1997| |Ducasse fc 



2003 



umir 



2008 ) and has been invoked when studying droplet collision rates in warm clouds ( Shaw 
Falkovich et q/.|[2002 ) and plankton encounter rates in the ocean ( Schmitt & Seu 



been observed in several experiments ( Aliseda et a/.|2002 


Wood et a/.|2005 


Salazar et al. 


2008 


Saw et a/.|2008 Monchaux et a/.|2010 


). Here we show direct experimental evidence 



that regions of high particle concentrations are also regions of high strain. 



2. Experimental method & observation 

We conducted three-dimensional Lagrangian Particle Tracking measurements (Ouel- 
|lette et aT]|2006| |Xu 2008) and obtained the trajectories of solid particles with different 
densities in a turbulent water flow. Below we give a brief summary of the experiment 
(see also Gibert et aL| ( |20f0 )). The turbulent flow was generated by two counter-rotating 
baffled discs of 28 cm in diameter. The turbulence chamber, shaped as a hexagonal cylin- 
der, measured 40 cm along the axis of the propellers and 38 cm in the cross-section, both 
in height (vertically) and width (horizontally). The axis of rotation of the discs was in 
the horizontal direction. The strong sweeping of the fluid near the bottom surface re- 
suspended the heavy particles with the effect that the particle concentration remained 
sufficient for the measurement of multi-particle statistics. A spherically shaped measure- 
ment volume with a diameter of 22 mm was located at the center of the water chamber, 
where the mean flow was small compared to the velocity fluctuations. The spatial resolu- 
tion of 45 /im per pixel and the temporal resolution of 31 frames per Kolmogorov time 
guaranteed accurate measurements of particle velocity and acceleration ( Xu et aL|2007 ). 
As shown in the Table[l]we used three different types of particles: polystyrene, glass, and 
stainless steel with densities between 1 ;$ Pp/pf ~ 8. For the Taylor microscale Reynolds 
number of the turbulent flow (R\ — 442), the particles had the size of the Kolmogorov 
scale (?7 = 74 /im). The resulting Stokes numbers, a relative measure of the particle iner- 
tia, were between 0.08 and 0.5 (see Table [l}. Other relevant parameters of the turbulent 
flow are the fluctuating velocity u' — 0.15 m/s, the integral length-scale L = 87 mm 
and the Kolmogorov time-scale t v = 5.2 ms. We measured 3 x 10 7 data points for the 
polystyrene particles, 8 x 10 7 for the glass particles and 3 x 10 6 for the steel particles. The 
fluctuations in Fig. [T] can be attributed to the statistical uncertainty. 
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Figure 1. Mixed Eulerian velocity-acceleration structure function {8 r v ■ S r a) normalized by the 
expected values for the fluid particles — 2e. The values of the different plateaus represented by 
the dot-dashed lines are : blue 1, red 1.15, black 1.30. 



For these heavy particles, we observed that the magnitude of the mixed velocity- 
acceleration Eulerian structure function (S r ~v ■ 5 r a. p } increased with particle density, where 
5 r v and <5 r a p are the relative velocity and the relative acceleration between two particles 
separated at distance r (e.g. 5 r v = v(x + r,t) — v(x,t)). For fluid tracers this mixed 
velocity-acceleration structure function in the inertial range equals a constant: (5 r u ■ 
5 r a) = — 2e, where u and a are fluid velocity and acceleration, respectively, and e is 
energy dissipation rate per unit mass (Mann et al. 1999 Falkovich et al. 2001; Pumir 



et al. 2001 Hill 2006). We also found this in our experiment, where the polystyrene 



particles behaved as tracers (Gibert et al. 2010). As shown in Figure [T] for inertial 



particles, this mixed velocity-acceleration structure functions in the inertial range were 
still constant, but the values increased significantly (up to 30%) when the particle Stokes 
number increased from 0.08 to 0.44. 

As reported earlier ( Gibert et a/.|2010 1 the relative velocities in this range changed by 
less than 10%. In addition, the RMS acceleration changed with particle density from 6.9, 
6.6 and 8.5 m/s 2 for tracer, glass and steel particles, respectively. Therefore, the observed 
change in (<5 r v • i5 r a p ) cannot be simply explained by changes in the magnitudes of <5 r v 
and <5 r a p , but must be associated with the alignment between the two vectors. Indeed, 
we observed a change of the cosines 



(La r 



|<5 r v||<J r a p | 



(2.1) 



of the angle between these two vectors with inertia of the particles. Considering the tracer 
particles (polystyrene), its average value (cos 8) monotonically increases from —0.052 to 
—0.01 when increasing the observation scale r (12 < r/rj < 120). For glass and steel 
particles, the average (cos#) increased in absolute value by 25% and 9%, respectively. 



3. Alignment of relative acceleration and relative velocity 

The effect of particle inertia on cos 6 is best demonstrated by its probability density 
function (PDF). Figure [2^l shows the PDFs of cos 9 corresponding to 5 r v and 8 r & p at 
separation r = (28 ± 7)r/ for the three different types of particles (similar behaviour 
is observed for other separations in the inertial range). All PDFs are skewed to the 
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Figure 2. (a) Probability density function of the cosine of the angle 9 between <5 r v and 8 r a p for 
the three types of particles for a separation r/77 = 28 ± 7. (b) Same PDFs as in (a) normalized 
by the PDF of the tracer particles, (b) shows the relative effect of the particle inertia on the 
alignment/anti-alignment of <5,-v with S r a p . 



side of cos 9 < 0, indicating that the two vectors are preferentially oriented in opposite 
directions, which gives the negative values of (S r v ■ S r a p ), as shown in Figure [T] 

A very interesting observation from Figure [2^, is the increase of the probability of 
I cos 6 1 ~ 1 with particle inertia. This can be seen even better in Figure [2Jd where the 
PDFs of the two types of inertial particles were normalized by the PDF of the fluid tracers. 
Clearly, the particle inertia enhanced the collinearity, either in the same direction or in 
the opposite direction, between the relative velocity and the relative acceleration vectors. 
We note that the PDFs are slightly skewed toward the anti-alignment side (cos 9 < 0). 

How can we understand this observation? For this range of Stokes numbers, as shown 



numerically and experimentally (Bee et al. 2006 Salazar & Collins 2011 Gibert et al. 



2010), the difference between particle and fluid velocity and acceleration at the particle 



position is very small. The same is true for the dynamical variables 5 r a p and 5 r v. Thus, 
if the inertial particles were homogeneously distributed in space, like fluid tracers, we 
should not observe any differences between the PDFs of cos for fluid tracers and inertial 
particles. This is not what we observed. Therefore our measurements strongly indicate 
that inertial particles preferentially explore certain regions of the flow, leading to the 
observed PDFs in Figure [2j 



4. Implication to preferential sampling 

To further elucidate which features of the turbulence are responsible for the observed 
behaviour, we investigate ^ r u and 5 r a of fluid tracers for simple flow fields. Let us consider 
a simple model, namely a steady, incompressible linear velocity field Ui defined by : 



Ui = MijXj 



(4.1) 



where tt, is the fluid velocity at position x i: My is the velocity gradient tensor. In our 
experiments, where the separation r are in the inertial range, we could treat My as an 
effective velocity gradient at scale r. It has been sho wn that the effective velocity gradien t 
at scale r evolves with a time scale to = {r 2 /e) 1 / 3 



Xu et al. 


2011 


Pumir et al. 


2011 



which is much larger than the particle relaxation time r p for r in the inertial range, as it 
can be quantified by noticing that the ratio r p /t = St(r]/r) 2 / 3 -C 1. Therefore, to study 
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Figure 3. PDF of cos 6 as defined in Eq. |2.f | The relative velocity and acceleration are obtained 
from two fluid tracers randomly distributed in a steady, purely straining linear velocity field. 
The different symbols corresponds to different ratios of A2/A1. The dashed-line corresponds to 
Eq.[0 



the dynamics of particle relative velocity and relative acceleration, the flow field may be 
treated as frozen. 

For a steady flow, the fluid acceleration field <2j is determined by My: dj = Uj?^ = 
MikMkjXj. The expressions for two-particle relative velocity and acceleration, together 
with the cosine of the angle between them, cos 9, can then be obtained in closed form. For 



a velocity field with pure rotation, M^, 



-Mji, it can be easily shown that 5 r a-5 r u = 0, 



i.e. the fluid acceleration is always normal to its velocity. The PDF of cos 9 is a ^-function 
at cos 6? = 0. This is not consistent with our observations. 

For a purely straining velocity field, the three eigenvalues of My are all real and can 
be arranged as Ai ^ A2 ^ A3. Incompressibility requires that Ai + A2 + A3 = 0. Then 
it is straightforward to show that the alignment between relative velocity and relative 
acceleration vectors depends only on the orientation of the separation vector r between 
the two particles and the ratio of the eigenvalues, e.g., A2/A1. The PDF of cos (9 at a given 
A2/A1 can then be obtained by assuming that r distributes uniformly in all directions 
within this purely straining flow. When A2 =0, i.e. the two-dimensional case, the PDF 
is: 

PDF(cos 9) = — 1 . (4.2) 

tt(1 — cos^ 9) L ' Z 

For more general cases, A2/A1 7^ 0, we obtained the PDFs numerically by choosing 
randomly 10 7 separation vectors r and computing cos 9 (Eq. 2.1| in the linear velocity 
field Ui (Eq. |4.1[ ). Please note, numerical analysis of various test cases showed that the 
shape of the PDF was insensitive to specifics of the particle distribution, ft is well known 
that in turbulent flows the average of the eigenvalue ratio of the velocity gradient is 
positive, with a value of (A2/A1) « 0.2 ( Betchov|1956 Ashurst et aL|1987 Kholmyansky 
et aL||2001 Liithi et aT|2005| |Berg et al. [2006 1 . For the effective velocity gradient in the 



inertial range, this ratio decreases monotonically from this positive value to when the 
scale r increases from rj to the integral scale ( |Pumir et aL|20lT| ) . We therefore performed 
simulations with the range of A2/ Ai around these values. Figure [3] shows the PDFs of cos 9 
from the simulation results with ^ A2/A1 ^ 0.3, together with the analytical solution 
Eq. 4.2 The PDFs for A 2 /Ai > in Fig. [3] are remarkably similar to the "relative PDFs" 
shown in Figure[2b, namely, the relative accelerations and the relative velocities prefer to 
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be collinear, with slight skewness to the anti-alignment side. This observation, together 
with the fact that inertial particles have, to leading order, the velocities and accelerations 
of the fluid at the particle positions ( Bee et aTf 2006; Sa lazar fc Collins|2011 Gibert et al. 



2010), suggests that inertial particles sample preferentially strain-dominated regions in 



turbulence. Our approach should be compared with results from numerical simulations 
of simplified equations, for which the phenomenon of "preferential sampling of strain- 
dominated regions" has been observed (Chun et al. 2005 Collins & Keswani 2004). 
There the velocity fields around the virtual inertial point particles was investigated 



It 



will be interesting to see how well the PDFs of cos 9 agree with the experimental results 
presented here. Moreover, studying the PDF of cos 9 for fluid tracers might reveal new 
information about the structures of the turbulent flow itself, e.g. by conditioning of the 
PDF on local strain rates and vorticities. 



5. Summary 

In summary, we observed a significant increase of the mixed velocity-acceleration struc- 
ture function (<5 r v • <5 r a p ) when increasing the inertia of heavy particles of the size of the 
Kolmogorov-scale. We found that this increase could not be attributed to the change in 
the magnitudes of the velocity and/or acceleration increments (\5 r v\ and/or |£ r a p |). A 
significant contribution to this effect is from the alignment or anti-alignment between the 
relative velocity and relative acceleration. Using a simple model, i.e. a linear velocity field, 
we show that in the inertial range heavy particles preferentially explore strain-dominated 
regions of a turbulent flow. 

The tendency of heavy particles to preferentially explore strain-dominated regions may 
also explain our recent observations (Xu & Bodenschatz 2008) that the time-averaged 
radial distribution function of particle positions increases when the center of the von 
Karman flow is approached. In that case, the long-time averaged mean flow imposes 
deterministic strain at the center of the apparatus. With our new findings reported here, 
one would expect that the long-time average would lead to a measurable preferential 
concentration towards the center of the apparatus. 

Our analysis does not rely on the commonly used and hard to measure long time 
average of the relative particle positions. Instead it only considers dynamical properties 
and thus opens a novel way to study particle turbulence interaction. For higher Stokes 
number, the filtering induced by the particle response time ( Bee et al| |2006; Salazar & 
CollinspOll ) will have an additional effect on the analysis presented above. 
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